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We consider the parity problem in one-dimensional, binary, circular cellular automata: if the initial 
configuration contains an odd number of Is, the lattice should converge to all Is; otherwise, it should 
converge to all Os. It is easy to see that the problem is ill-defined for even-sized lattices (which, by 
definition, would never be able to converge to 1). We then consider only odd lattices. 

We are interested in determining the minimal neighbourhood that allows the problem to be solv- 
able for any initial configuration. On the one hand, we show that radius 2 is not sufficient, proving 
that there exists no radius 2 rule that can possibly solve the parity problem from arbitrary initial 
configurations. On the other hand, we design a radius 4 rule that converges correctly for any initial 
configuration and we formally prove its correctness. Whether or not there exists a radius 3 rule that 
solves the parity problem remains an open problem. 

1 Introduction 

Understanding the nature of computations within cellular automata remains an elusive problem. In fact, 
in spite of their long-proclaimed ability to perform computations, very little is still known as to how we 
should design the local state transitions towards achieving a given global behaviour. As examples are 
designed or found by search, it is inevitable to try to understand their underlying programming language; 
but the truth is, to this date, every attempt along these lines has fallen into the strenuous effort of trying 
to tame local state patterns towards the global state target, or trying to make sense of the latter in terms 
of the former 0. 

On the other hand, studying how to employ local actions to achieve desirable global behaviours 
is of utmost importance and extensively investigated in many other evolving systems (e.g., distributed 
systems, mobile robots, population protocols). In such systems, in fact, understanding the limitations and 
the power of local interactions to solve global computations has immediate implications on the design 
of efficient and scalable solutions (e.g., see HI El [03 El). Cellular automata (CAs) are the simplest 
possible evolving systems, and understanding the impact that neighbourhood size has on computability 
could have consequences for more complex systems based on local interactions. 

This paper is aligned with these efforts. Here, we concentrate on the one-dimensional parity problem, 
which has essentially the objective of figuring out the parity of an arbitrary binary string, by means of 
a one-dimensional, binary cellular automaton fl4l . The parity problem is a well-known benchmark 
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task in various areas of computer science, typically camouflaged under the XOR operation on a binary 
input, as in artificial neural networks [7], but it also lends itself to the context of cellular automata, as 
a typical case of a global problem that has to be solved by purely local processing. The problem is 
formulated under periodic boundary conditions and arbitrary finite lattice size, so that, if the parity of the 
global configuration is odd, the lattice is supposed to lead to an homogeneous configuration with only 
Is; otherwise, it should converge to all Os [9]. 

The notion of parity has appeared quite often in the CA literature, even if implicitly, as it bears 
relevance to the related notion of additivity of CA rules |3l[T5]]. However, the parity problem per se 
has not been extensively studied, particularly in comparison with the well-known benchmark CA task of 
density classification, where the aim is to determine the most frequent bit in the initial configuration of 
an odd-sized lattice, also by reaching an homogeneous configuration. The density classification problem, 
in fact, has been extensively investigated and fully understood in odd-sized lattices. In particular, it has 
been shown that there exists no single rule able to solve the problem for any arbitrary initial configuration. 
Combinations of rules have been devised, however, as well as probabilistic solutions to the problem (e.g., 

see da annua). 

An advantage in favour of the parity problem is that, from the perspective of automata theory, it is 
simpler than its kin, insofar as the notion of parity can be handled by finite automata, whereas the ability 
to compare arbitrarily variable quantities (which is inherent to density classification) requires at least a 
pushdown automaton JSJ. In fact, the increased simplicity of the parity problem is reflected in the fact 
that it is easier to find good rules for it, by searching, than for density classification ifTol . Therefore, there 
are strong reasons for considering the parity problem generally more tractable and amenable to analysis, 
which makes it a serious candidate for case studies that might help the understanding of the nature of 
computation in CAs in general. 

The parity problem is ill-defined for even-sized lattices (by definition, an all 1 configuration con- 
verges to an all configuration making it impossible for any rule to converge to 1). Modifying the 
definition of the problem to allow the target homogeneous configuration to be achieved only once, and 
not as a fixed point, the problem becomes solvable also for even lattices. In fact, by relying on this 
variation, it can be perfectly solved by a carefully engineered sequence of rule applications, quite sur- 
prisingly, of elementary cellular automata ifTTTl . However, if we do not want to change the definition of 
the problem, it is then necessary to restrict the study to odd-sized lattices. We then say that a CA rule is 
perfect if it solves the parity problem for arbitrary initial odd-sized configurations. 

Unlike the density classification problem, we show that the parity problem can, indeed, be solved 
by a single rule. Besides being interested in its general solvability, we are also interested in determining 
the minimal neighbourhood that allows the construction of a perfect rule. With this goal in mind, we 
first prove that radius 2 is not sufficient for a perfect rule to exist. We first identify several constraints to 
which such a perfect rule is subject and we show that no rule is feasible with all of these constraints. We 
then show that the problem becomes solvable when CAs have radius 4: our proof is constructive as we 
design a perfect rule and we prove its correctness. We leave open the case of radius 3, for which there 
is strong empirical evidence that no perfect rule exists, but that there might be radius 3 rules that would 
solve the problem for prime- sized lattices iPToll . 

2 Notation and Basic Facts 

We consider one-dimensional, binary cellular automata (CAs) on finite lattices with periodic boundary 
conditions. Let / : {0, 1} +1 — > {0,1} denote the local rule of a CA with radius r. The global dy- 
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namics of a one-dimensional cellular automaton composed of n cells and radius r is then defined by the 
global rule (or transition function): F : {0, 1}" -> {0, 1}" s.t. VX £ {0, 1}", Vi 6 {0, . . . ,n- 1},F(X),- = 
f(x[- r ...,Xi,... ,Xi+ r ), where all operations on indices are modulo n. 

A fixed point P G {0, 1}" of a circular CA with global transition rule F is a configuration P such that 
F(P) = P. 

We say that a cellular automaton converges to a configuration P from configuration X° if P is a 
fixed point and if for some finite n, F"(X°) = P where F" is the « tn iteration of F. We are particularly 
interested in the homogeneous configurations as fixed points and will refer to these as the 0-configuration 
and the 1 -configuration. 

We say that a local rule solves the parity problem if, starting from an arbitrary initial configuration, 
on an arbitrarily sized lattice, the cellular automaton converges to the 0-configuration, if and only if the 
initial configuration contains an even number of Is, and converges to the 1 -configuration otherwise. 

Since a rule solving the parity problem must converge to the homogeneous configurations, we have 
our first two properties, of perfect rules. 

Property 1 If f solves the parity problem, then /(0, • • • , 0) =0 and /(l , ■ • • , 1) = 1. 
It is immediately obvious that, by definition, no solution exists for even-sized lattices. 
Theorem 1 Consider circular CAs with radius r and even size n. There exists no rule that works cor- 
rectly from any initial configuration. 

PROOF Trivially /(l . . . 1) =0 otherwise the configuration with all Is would incorrectly converge to 1. 
Since (11111. ..Ill) is not a fixed point, it follows that no initial configuration can ever converge to the 
1 -configuration, g 

For this reason, from now on, we consider only odd-sized lattices and we call a rule perfect if it solves 
the parity problem for any odd-sized lattice, starting from any initial configuration. 

We now recall the definition of de Bruijn graphs, which are useful tools for representing CA rules and 
which will be used in the subsequent sections. The de Bruijn graph of a local rule of radius r is a directed 
graph on 2 nodes, one for each value in the set {0, l} 2r . There is an edge from node xq . . .X2,—i to node 
y\ ...yi r if Xj = yt for all i from 1 to 2r — 1. These edges are labelled with the value of the local function 
at (xq, . . . ,x 2r -i,y2r), f(xo, ■ • • ,x 2r -i ,y2r) = f(*o,y\ , • • • ,yir)- Note that the shape of the de Bruijn graph 
for a local rule of a given family (i.e., those with the same neighbourhood and states) is fixed, only the 
edge labels change. For example, Figure [T] shows the de Bruijn graph for a radius 1 rule. 




Figure 1: De Bruijn graph for the local parity rule (150). 
The following is a necessary condition for parity preservation. 
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Property 2 A rule preserves the parity of a configuration if active transitions always come in pairs. 
That is, given a local rule f of radius r, and any configuration (xq,... the number of times that 

f{xi-„ ...,Xi,... ,x i+r ) / Xi is even. 

It is also very simple to see that no solution exists for elementary circular CAs (i.e., with radius 1). 
Theorem 2 There exists no perfect rule for elementary CAs. 

PROOF From Property [I] for any perfect rule, we must have /(000) = and /(111) = 1. Now con- 
sider a configuration containing a single 1 . In order to both maintain parity and move towards conver- 
gence, we must have /(100) = /(010) = /(001) = 1. Similarly, from the singleton 0, we must have 
/(110) = /(101) = /(Oil) = 0. So the only possible perfect rule is the local parity check (rule 150). 
However, it is easy to see that such a rule does not converge from several initial configurations, for ex- 
ample from (.. .0001000. . .). | 



3 Impossibility with Radius 2 

In this section we show that with radius 2 it is impossible to construct a perfect parity rule. 

Our aim is to show, first, several necessary transitions for a perfect rule, and second, the existence 
of a limited set of feasible pre-images for the two final homogeneous configurations. Each possible pair 
of feasible pre-images further induces necessary transitions for a perfect rule, significantly reducing the 
space of possible perfect rules. We conclude the proof by running exhaustive tests to verify that in this 
set there exists no perfect rule having the necessary transitions. We begin with a series of lemmas that 
force certain transitions to or 1 . 

Consider the de Bruijn graph for radius 2 rules. A pre-image of the final 0-configuration corresponds 
to a cycle of odd size and even parity. A pre-image of the final 1 -configuration is a cycle of odd size and 
odd parity. Let 38$ be the subgraph containing only the edges corresponding to transitions to and 8§\ 
the subgraph containing the edges corresponding to transitions to 1 . 

Lemma 1 In a perfect rule, three or five of the following must transition to 1: 
(10000), (01000), (00100), (00010), (00001). 

PROOF A configuration consisting of a single 1 must eventually converge to all Is, hence the number of 
Is in the configuration must increase. Furthermore, in order to maintain parity, it must increase to an odd 
number. The five configurations above are the only ones occurring at the local level that are not all 0, and 
therefore 3 or 5 of them must go to 1 . ■ 

Similarly, 

Lemma 2 In a perfect parity rule, three or five of the following must transition to 0: 
(01 111), (10111), (11011), (11101), (11 110). 

Lemma 3 Neither 83$ nor 88 \ can contain a cycle of even length and odd parity. 

PROOF A cycle of even length in either 83® or 88\ will become a sequence having even parity at the next 
iteration since it will be either all 0s or an even number of Is, so this cycle itself will have changed parity. 
Assume that the de Bruijn graph of a rule F admits such a cycle, let C be such a cycle and let P be a cycle 
of odd length passing through a node of C. For the rule to be perfect, F (P) must have the same parity as 
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P. Now consider a new cycle P' formed from P by adding the cycle C where P passes through it. Since 
C has even length, P' has odd length. Since C has odd parity, P and P' have different parity. However, 
since F(C) has even parity, F(P') = F(P), hence the parity of P' has changed and F cannot be perfect, g 

Lemma 4 In a perfect parity rule either: i) /(10101) = 1 and /(OIOIO) = 0, or ii) /(10101) = and 
/(01010) = 1. 

PROOF This is a direct consequence of Lemma [3] since /( 10101) = /(01010) would imply the existence 
of the even cycle with odd parity (1010,0101, 1010,0101, 1010,0101) either in m or in SS\. a 

Lemma 5 In a perfect parity rule, it is impossible to have four or more consecutive 0s in a pre-image of 
the 0-configuration. 

PROOF Let the pattern (0000) be present in a pre-image P of the 0-configuration. Then the follow- 
ing neighbourhood configurations must also be present and must be transitioning to 0: (10000) and 
(00001). By Lemma [T] we must then have the following configurations transitioning to 1, so they 
may not occur in P: (01000), (00100), (00010). Hence our group of four 0s must be both pre- 
ceded and followed by at least two Is, thus entailing that we have the necessary transition set 5 = 
{/(11000),/(00011),/(10000),/(00001)} -> 0. Consider now an initial configuration containing a 
single 1 surrounded by 0s. From S, we have that 0001000 can only grow to 001 1 100, but, again from S, 
we have that, from 0011 100, no growth is possible anymore, which is a contradiction. | 

Analogously, we have: 

Lemma 6 In a perfect parity rule, it is impossible to have four or more consecutive Is in a pre-image of 
the I -configuration. 

From Lemma [3} any feasible pre-image of the 0-configuration (resp. 1 -configuration) is either a 
simple odd cycle c with even (resp. odd) parity, or the composition of cycles not containing any even 
cycle of odd parity. 

So, to identify feasible pre-images for final configurations for lattice size n in the de Bruijn graph, 
we have to find at least one cycle of size n to be labeled and one to be labeled 1, having the property 
that they do not include: 

(i) the self-loops (which are forbidden by Lemmas [5] and [6] 

(ii) the 2-cycle (0101,1010) (which is forbidden by Lemma[4j); and 
(Hi) an even cycle with odd parity (Lemma[3]). 

By inspecting all cycles of size 5, we obtain that: 

Lemma 7 In a perfect rule at least one of these three cycles in the de Bruijn graph must transition to I: 
B\ = (0011,0111,1110,1100,1001) (corresponding to configuration: 00111) 
B\ = (0000,0001,0010,0100, 1000) (corresponding to configuration: 00001) 
B\ = (0101,1011,0110,1101,1010) (corresponding to configuration: 01011) 
and one of these must transition to 0: 

= (0001,0011,0110, 1100, 1000) (corresponding to configuration: 00011) 

= (0111,1111,1110,1101,1011) (corresponding to configuration: 01111) 

= (0010,0101, 1010,0100, 1001) (corresponding to configuration: 00101) 
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PROOF B 5 { , B\ and B\ (resp. , W% and W 3 5 ) are the only cycles corresponding to feasible pre-images 
for the 1 -configuration (resp. O-configuration) for lattices of size 5, which do not violate Lemmas |4j|5] 
and [6] ■ 

Consider, now, lattices of size 7. All cycles of length 7 have been enumerated and the only cycles 
that do not contradict Lemmas|3]|4][5] and[6]and correspond to feasible pre-images of the 1 -configuration 
are: 

B\ = (0000, 0001 , 001 1 , 01 1 1 , 1 1 10, 1 100, 1000) (configuration: 00001 1 1) 
B\ = (0001,0011,0110,1101,1010,0100,1000) (configuration: 0001101) 
B] = (0001,0010,0101,1011,0110,1100,1000) (configuration: 0001011) 
B\ = (1001,0011,0110,1100,1001,0010,0100) (configuration: 1001100) 

Analogously, the only cycles which do not violate Lemmas [3] |4j [5J and [6] and correspond to feasible 
pre-images of the 0-configuration are: 

Wj = (0001,0011,0111,1111,1110,1100,1000) (configuration: 0001111) 
W 7 = (0010,0101,1011,0111,1110,1100,1001) (configuration: 0010111) 
= (0011,0111,1110,1101,1010,0100,1001) (configuration: 0011101) 
Wj = (0110,1100,1001,0011,0110,1101,1011) (configuration: 0110011) 

From simple observation, we can rule out some of these cycles and combinations of cycles. 

Lemma 8 A perfect rule of radius 2 cannot have as a pre-image of the 0-configuration. 

PROOF Cycle shares at least one transition in common with each of the possible pre-images of the 
1 -configuration of size 7. For example, W^, B\ and B\ all share the edge (0001,0011) in the de Bruijn 
graph. Cycles W? and 5^ share (0010,0101), and W? shares (0011,0110) withS^. a 

Lemma 9 A perfect rule of radius 2 cannot have B\ as a pre-image of the 1 -configuration. 

PROOF First, if B\ is a pre-image of the 1 -configuration, then W 2 5 is a pre-image of the 0-configuration 
of size 5 since B\ and W 3 5 share (1010,0100). Cycle B] also has transitions in common with Wj, 
Wj 7 , Wj. It has no common transitions with W^J, however, W 2 7 and W 2 5 together form a cycle, namely, 
(01 11, 11 10, 1101, 1011), in violation of Lemma|3] a 

Similarly, we can show, 

Lemma 10 A perfect rule of radius 2 cannot have B^ as a pre-image of the 1 -configuration. 
In fact, we can now restrict to very well defined possible cases. 

Lemma 11 A perfect rule of radius 2 must have W 2 as a pre-image of the 0-configuration, B\ as a 
pre-image of the 1 -configuration and either 

• B 7 { as a pre-image of the 1 -configuration and Wj as a pre-image of the 0-configuration, or 

• B 7 4 as a pre-image of the 1 -configuration and Wj 7 as a pre-image of the 0-configuration. 
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PROOF From the lemmas above, we know that the only possible pre-images of the 1 -configuration of size 
7 are B\ and B\. Cycle B\ has transitions in common with all possible pre-images of the O-configuration 
except Wj. Of the 5-cycle pre-images of the 1 -configuration, only B 5 2 is compatible with Wj. Of the 
possible pre-images of the O-configuration of size 5, neither W 2 5 nor poses any conflict, however, W 3 5 
makes it impossible to have any cycles of length 3 going to the 1 -configuration, so that all configurations 
having a period of size 3, (i.e. configurations of the form (001001001 • • • 001) will fail to converge). The 
proof is analogous beginning with cycle B\. | 




Figure 2: Possible perfect rule for radius 2, with B\ and Wj. 
We can finally conclude: 

Theorem 3 There is no perfect parity rule of radius 2. 

PROOF From the lemmas in this section it follows that for a perfect parity rule of radius 2, we must have 
one of the de Bruijn graphs shown in Figure [2] or Figure [3] Note that, in these graphs, depending on 
the choice of xi, we can place further restrictions on the remaining unknown edges (again, due to the 
lemmas in this section). Consider Figure |2j if x\ = 1, then in order to avoid creating cycles with an odd 
number of transitions, we must have a\^c\,b\j= d\ from cycles (0001,0010,0101, 1010,0100, 1000) 
and (0101,1011,0111,1110,1101,1010). Whenjci = by contrast, we must have a x = c x ,b\=d\. The 
constraints on ai, bi, ci, and d2 resulting from the choice of xj are the same. Testing has shown that all 
16 of the resulting rules will fail for some initial configurations. | 

A final note concerns lattices of prime size. It has been conjectured, in the case of radius 3, that there 
may exist rules with the desired behaviour on arbitrary lattices of prime size Ifl6l . While the case of 
radius 3 is still open, we can show that the impossibility result for radius 2 holds even if we restrict the 
discussion to prime-sized lattices. 

Theorem 4 There is no radius 2 rule the always solves the parity problem even restricting to lattices of 
prime size. 
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Figure 3: Possible perfect rule for radius 2, with B\ and Wj. 



PROOF Restricting to prime size lattices, we can no longer use Lemma 1 1 This introduces only a few 



possible extra cases using B 3 or W 3 . Consider rules containing By As before, we can eliminate W±, 
B\ and B\ as pre-images of the 0- and 1- configurations. In addition, B\ conflicts with IV 2 7 and W 3 on 
the edge from 0101 to 1011 and with Wj on the edge from 0110 to 1101. Hence, we must have Wj as 
the 7-cycle pre-image of the 0-configuration. Now, B\ conflicts with Wj, so we are left with B\ as the 
7-cycle pre-image of the 1 -configuration. Since B\ conflicts with W 3 on the edge from 0100 to 1001, 
we are left with W 2 5 as the 5-cycle pre-image of the 0-configuration. These results are illustrated in the 
graph of Figure[4] Proposition[2]dictates that the edges labeled c?3 must be the same, as will be the edges 
labeled b^. Now consider rules containing W 3 5 . Similar analysis shows that must have B%, Wj and B\, as 
illustrated in Figure [5] Testing has shown that all 8 of these rules will also fail for some initial configu- 
rations of prime length. | 



4 A Perfect Rule with Radius 4 

We now describe the construction of a rule with radius 4 having the desired properties: parity preservation 
and convergence. The intention is to first give the reader an intuitive understanding of how the rule works 
and how it was developed. Formal proofs will follow in the subsequent section. 

4.1 Rule BFO 



The most compact representation of rule BFO that we propose for solving the parity problem in any 
lattice of odd size is given in Figure [6] and corresponds to rule number: 

127660 1 9579927887748828308783632 1 25 1 3720894862957 1434 19940439400267 1 69599 1 869267727- 
072917454377539194754200976283425175983876539715064584172642413634846720 
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Figure 4: Additional possible perfect rules for radius 2 on prime lattices with Wf. 




Figure 5: Additional possible perfect rules for radius 2 on prime lattices with B\. 

in Wolfram's lexicographic ordering scheme. The figure shows all active transitions (i.e., transitions that 
change the current state). However, it is often easier to explain why and how the rule works using a less 
compact form, where pairs of rules can be made explicit. This form of the rule is given in Figure |7J we 
will be referring to this representation in the remainder of this section. 

We now describe the intended behaviour of the rule before proving its correctness. Consider an initial 
configuration X° as being formed by blocks b\ of consecutive Is separated by blocks w\ of consecutive 
Os: X° = (b^w® ,^2)^2) " " ibk> w k " ')• The idea of our construction is to have a block bi of Is propagate 
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* 1 1 100*** 
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*00100*** 
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00100**** 
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**010100* 


1 


1 1 101**** 





*0101*0** 





**01 10*** 





***110110 





***0110** 





****1101* 






Figure 6: Minimised rule BFO (the asterisk * refers to any value). 
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T 3 ,T 4 
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T 5 ,T 6 


Annihilation of 1 1 
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Local shift 


T 6 : 


**0110*** 







?9,7lo 


Leftward growth of 0-blocks 




*010100** 









Start of 0-growth 


r 9 : 


***11101* 









Local adjustment 


2io: 


111010*** 









?ii: 


1110111** 









Tn: 


**1110110 










Figure 7: Active rule transitions (left) and behaviour of combinations of rules (right). 

to the right, two cells per iteration, until a stopping condition or convergence has been reached. Such 
propagation might result in merging the block with the next bi+\ (if the corresponding w; is of even size). 
When the merger does not occur (because |w;| is odd or due to some other condition), there will be a 
propagation of 0s to the left, led by a block of the form (01). Such counter-propagation might result in 
the total annihilation of the block of Is. Otherwise, it will result in the creation of a single 1 surrounded 
by 0s, which will start propagating to the right again. We will show that such behaviour reduces the 
number of blocks, eventually converging to an homogenous configuration. 

We now describe some properties of the rule that can easily be derived by construction and that give 
an intuition for the reasons for the behaviour described above. Note that, by construction, the rule's 
transitions always occur in pairs; in other words, whenever a transition occurs in a cell, another transition 
occurs in its neighbourhood. It is useful to describe the behaviour of each pair and we will also use these 
pairs to prove that parity is being preserved. 

— Rightward growth of 1-blocks. A singleton 1 grows to the right, by two Is at each step, if it is 
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preceded and followed by at least two 0s, as prescribed by transitions Tt, and T4. A block of three or more 
Is grows to the right if it is followed (on the right) by at least two 0s. This behaviour is created by the 
pair of transitions T\ and T2. 

— Annihilation of pairs of Is. As a consequence of transitions T$ and T$, an isolated pair of Is is always 
eliminated. 

— Leftward growth of 0-blocks. A (01) block moves to the left, leading a growing block of 0s (at a 
growth rate of two 0s per step) if there are at least three Is to its left and one of the following: (i) at least 
three Is to the right of the (the growth is obtained by the pair of transitions, (Tg,Tu)', or (ii) at least one 
to its right (due to Tg and T\o). Note that the pair (Tg, T\\) starts the growth of a 0-block, while the pair 
(Tg , Tio) continues the growth as far as possible. 

— Local Shift. A (101) block is transformed into (110) if there are a on its left and at least two 0s on 
its right (combination of transitions T-] and r§). 

— Local Adjustment. Finally, if a (01 10) block is preceded by at least three Is that is, (...1110110...) 
occurs, in order to avoid parity errors due to the annihilation of the pair of Is, we force the creation 
of a solid block of 0s to the right of the existing block of Is with transition pair {T^,T\2), so that 
(...1110110...) becomes (...1000000...). 

Examples of the evolution of the rule are given in Figure [8] 




Figure 8: Evolution of rule BFO for even parity (left) and odd parity (right). A black cell corresponds to 
1, a white cell corresponds to 0. The initial configuration is at the top and time goes downward. 

4.2 Correctness 

In order to show that the rule we have constructed (or in fact any rule) performs a perfect parity check, 
we must prove that it preserves parity at every iteration, and that it always converges in finite time to an 
homogeneous configuration. We begin with the proof of parity conservation. 

4.2.1 Parity Preservation 

A rule preserves the parity of a configuration if active transitions always come in pairs. That is, given a lo- 
cal rule / of radius r, and any configuration (xq, • • • ,x«-i ), the number of times that /(x,_ r , . . . ,x,-, . . . ,x,- +r ) 7^ 
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Figure 9: Reduced transitional de Bruijn graph for rule BFO. 

Xi is even. To show that our rule does indeed have this property, we will use a modification of the de 
Bruijn graph. 

Given a configuration X = (jci, . . .x n ), one can determine the next iteration, F(X) by reading the edge 
labels as one traverses the graph from (jc,-_4 . . . x,- . . .jCi+3) to (jc,-_3 . . . x, . . .jc,-+4). Since we are considering 
circular configurations, the traversal of the de Bruijn graph will result in a closed loop. For our purposes, 
we are not interested in what the actual output is, only if an active transition has occurred. In keeping 
with that, we are, in fact, only interested in the parts of the graph where such transitions occur. Since a de 
Bruijn graph for a function of radius 4 can be quite unwieldy, we define a reduced transitional de Bruijn 
graph modifying the standard de Bruijn graph as follows. First, we label the edges with 7} (indicating 
one of our 12 active state transitions), or N, meaning no transition. Second, we draw only those parts of 
the graph connected to transitions, reducing the rest of the graph to a single node, denoted by an asterisk. 
Also, where there is no conflict, several nodes leading to or from the same transition are depicted as 
one, using the notation of the previous section; for example, there is an edge from node (*11100* *) to 
(11100* **), because of the various state transitions entailed by T\, given all possible values for the * 
symbols. Furthermore, we ensure that all nodes in the reduced transitional de Bruijn graph are distinct 
for all values of the * symbol. Finally, for any nodes occurring explicitly in the graph, all adjacent edges 
are represented, whether they correspond to a transition or not. It is easy to see that parity preservation 
can be detected from the reduced transitional de Bruijn graph for a given rule; more precisely: 

Lemma 12 A rule is parity preserving if and only if any cycle in its reduced transitional de Bruijn graph 
contains an even number of edges labeled with some transition 7}. 

PROOF Since the output of a circular CA is given by the edge labels of a cycle in its de Bruijn graph, we 
need only count the transitions to verify parity preservation. Furthermore, the reduced graph compresses 
only parts of the graphs where no transitions occur. Hence, if the are no cycles containing an odd number 
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of transitions in the reduced graph, there can be no configurations leading to an odd number of transitions 
and vice versa. | 

By inspecting the transitional de Bruijn graph for rule BFO, and by noticing that there are no cycles 
containing an odd number of transitions, we then have: 

Theorem 5 Rule BFO is parity preserving. 
4.2.2 Convergence 

We now turn our attention to the more challenging problem of proving that this rule will converge under 
any conditions. We can think of any CA configuration as an alternating sequence of blocks of Os and 
blocks of Is of varying lengths. We will show that BFO eventually converges by showing that its only 
fixed points are the homogeneous configurations and, furthermore, that any change in the configuration 
will lead, in a finite number of iterations, to a reduction in the overall number of blocks. Our first lemma 
shows that every non-homogeneous configuration is changing. 
Lemma 13 The only fixed points of rule BFO are the homogeneous configurations. 
PROOF A fixed-point configuration cannot contain pairs of Is, since rule pair (Ts,Ts) would apply. It 
cannot contain a block of three or more Is since (7i,72) would apply if it is followed by two or more 
Os, and a pair containing T 9 would apply if it is followed by only one 0. Therefore, a non-homogeneous 
fixed-point configuration could only contain isolated Is but the odd length of the configuration would 
imply that we must have at least two consecutive Os, hence the sub-configuration 0100 must occur. If 
this is preceded by a 0, (73,74) apply. If it is preceded by a 01, (Tj,T$) apply. | 

We now wish to show that every transition pair will eventually lead to a reduction in the total number 
of blocks. We begin with the transition pairs for which this is immediate. 

Lemma 14 Transition pairs (T$^), (T],T&) and quadruplet (T9, 7i2, T$,T^) lead to block reduction in a 
single step. 

Lemma 15 Transition pair (T 9 ,T\q) leads to block reduction in a finite number of steps. 

Proof Rule pair T 9 : /(* * *11101*) = and T U) : /(ll 1010 * **) = causes the leftward growth of 0- 
blocks. While a single application of this pair maintains the number of blocks, it leads (possibly through 
a repeated application of the pair) to an eventual block reduction through either an annihilation or the 
creation of a single 1 : 

01111010- •• ~> 01101000- •• by rules T 9 and T l0 
~> 00001000 • • • by rules T 5 and T 6 
0111010- •• ~> 0101000- •• by rules T 9 and T 10 
~> 01 10000 • • • by rules T 7 and T% 
~> 0000000 • • • by rules T 5 and T 6 



Notice that (T9, 7\o) leads to reducing the number of blocks by either two or four, depending on the parity 
of the block of Is on which it is acting. Also note that, even though it is possible for the leading block 
of Is to have shrunk from the left side, while the (T 9 ,T\o) pair is reducing it from the right, one of these 
two situations will still be reached, since, on the left, the Is can only be eliminated one at a time. ■ 
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Lemma 16 Transition pairs (T\ , T%), and (T%^T\) lead to either reduction or maintenance of the number 
of blocks. 

Proof Both pairs T\ : /(*11100***) = 1, T 2 : /(11100****) = 1, and T 3 : /(*00100***) = 1, 
T4 : / (00100 ****) = 1 are responsible for the rightward growth of 1-blocks. In fact, they grow a block 
of Is until it either merges with the next block or an isolated preceded by three or more Is is created. 
At this point, one of the following transition sets will apply: (T9, 7io) if the isolated is followed by 10, 
(T$,T\2, Ts,T(,) if it is followed by 1 10, and (Tg,Tu) if it is followed by three or more Is. We have already 
seen that the first two cases lead to block reduction, only the latter case can maintain block numbers. | 



Lemma 17 Transition pair (7g , 7\ 1 ) leads to reduction or maintenance of the number of blocks in a finite 
number of steps. 

Proof The pair T 9 :/(***11101*) =0, T n :/(*** 11 101 1) =0 is responsible for the start of 0-growth. 
This is the only rule pair that initially increases the number of blocks. It is the beginning of the growth 
of 0s from an isolated surround by three or more Is on either side. Once this growth has begun, it 
is continued by the (T9, T\q) pair until the number of blocks has been returned to its original size or is 
reduced by two. | 

Note that in the proof above, if the original number of blocks is maintained, it is because the tran- 
sitions have produced an isolated 1 with two or more 0s on either side which then begins to grow to 
the right either merging with the block on the right (and thus reducing the total number of blocks) or 
creating an isolated which then begins to grow left. What we wish to avoid is a CA which evolves to 
a periodic configuration of Is growing until only an isolated remains and then shrinking back to an 
isolated 1 which then regrows. We call this pattern of growing and shrinking the accordion effect. In the 
next lemma, we show that the accordion effect cannot occur on lattices of odd length. 




Figure 10: Accordion loop. 



Lemma 18 The accordion effect can only occur on lattices of even length. 



PROOF Figure 10 shows the transition pairs involved in a cycle of growing and shrinking. If we think of 
this loop as starting from the (Tg,Tn) node, a block of Is is shrinking until the (T$, T^) node is reached 
and then beginning with the (73,74) pair, a block of Is begins to grow. Let us assume that the entire 
lattice is perpetually in some stage of this cycle. We make several observations: 
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• The 1 -blocks start off having odd length. 

• The 0-blocks start off having even length. 

• In order for a block of Is to regrow, it must have even length at the end of the shrinking process. 

• In order for a 0-block to grow, it must have odd length at the end of the growth of Is. 

Since the entire configuration is in this process, it is made up of sub-configurations of two forms: 
a block of Is followed by a block of Os when the Is are growing; or bfilwi a block of Is separated by a 
01 from a block of 0s when the block of Is is shrinking. Consider a block of 0s, w\ followed by a block 
of Is, hi- If w\ has odd length, then it must not change lengths again before the growth of 0s restarts 
otherwise when the block of Is on its left, b\ grows, it will merge with bi. This means that (T9, Tn) must 
be applied on the left before (T$, T(,) can be applied on the right. Hence bi will shrink by 1 before the Is 
have finished shrinking on the right. Since we will need to have an even number of Is at that time, we 
must now have an odd number of Is. In other words, a block of 0s of odd length is always followed by a 
block of Is of odd length. Now assume that b\ has even length. In this case it must change sizes before 
the regrowth of Is is complete, so (T5, T&) must be applies on the right before (Tg, Tn) is applied on the 
left again. Hence W2 must already have even length. Taken together, we see that for the accordion effect 
to endure in perpetuity, the CA must have even length. | 




Figure 1 1 : The accordion effect on an even-sized lattice (left) and the successful resolution of an odd- 
sized lattice (right). 

Finally, taking these various lemmas together, we have the heart of our convergence proof. 

Lemma 19 From any non-homogeneous configuration, the total number of blocks decreases in finite 
time. 

PROOF From the previous lemmas, we see that only rule pair (Tg,Tn) increases the number of blocks 
and that within a finite number of steps, this increase is resolved. Our contention is that the accordion 
effect is the only way to maintain block numbers. Consider a {T^,T\\) pair occurring anywhere in our 
CA. If the leading block of Is is odd at the end of the execution of the (79, T\q) pairs, then the number 
of blocks had decreased by two. If it is even, the Is will begin to grow right. If the 0-block created by 
the execution of (Tg,Tn) is unchanged from the right, then it has even length and the regrowth of Is will 
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result in a merger of Is and a reduction in the number of blocks. Now the only way for the 0-block to 
have changed on the right is if the 1 -block on its right had shrunk due to the application of {T^,T\q) or 
(79,7ii). If the reduction in the 1 -block was not initiated by the (T$, 7n) pair, then real reduction in the 
number of blocks has occurred. If it was initiated by (Tg, T\\), then block reduction can only be prevented 
if another (79, T\\) pair is being executed to its right. Arguing in this way, we see that reduction in total 
block number can only be avoided if the CA is experiencing the accordion effect which can only occur 
in even-sized lattices, g 

From Theorem [5] and Lemmas [T3] and [19] we obtain: 

Theorem 6 Given a CA of odd length, rule BFO converges to all Is if the initial configuration has odd 
parity and to all Os if it has even parity. 

5 Concluding Remarks 

In this paper, we have established upper and lower bounds on the radius of rules that solve the parity 
problem by showing that there exists a rule of radius 4 which converges to all Is if the initial configuration 
is odd, and to all Os if it is even and, further, by proving that this problem is unsolvable by rules of radius 

2, even with the less strict condition of prime-sized lattices. The corresponding questions for radius 3 
remain open. However, we have developed tools in this paper that should be helpful in solving this latter 
problem as well. 

It is clear by now, how painstaking the task of designing a CA rule can be, let alone the formal proof 
of its correctness. To some extent, the process reminds us of similar programming efforts on simple, 
pre-modern computational models, such as Turing machines. And in this sense, we are still indeed at 
this point in history, when programming cellular automata. 

Since our main motivation for addressing the parity problem is not conscribed to it, one may ask how 
generalisable our experience herein could be to related problems, including the parity problem for radius 

3, as well as other related computational problems for CAs. It is tempting to think of the possibility of 
implementing a high-level programming approach that would automatically generate the state transitions 
of a CA rule, given the kinds of notions we have used, such as the growth of blocks of a given size in 
a given direction, the annihilation of blocks of given kind, etc. Even if this form of programming, so to 
speak, by patterns, does not solve the problem of designing a rule (the target algorithm), at least it would 
help its high-level conception, and its implementation in terms of the required state transitions. 

As a methodological note, it is worth mentioning that it was demanding in practice to resort to 
computational aids to complement the formal efforts. After all, the details involved in rule design are so 
many that it is quite easy to overlook some of them. This turned out to be essential in the present case for 
fine tuning our design in its origin. Such an interplay between formal and computational methods also 
came into play for devising the most compact representation of the BFO rule, as shown in the paper, for 
enumerating all cycles of length 7 in the de Bruijn graph of radius 2, and for the evaluation of the radius 
2 rules that had retained potential for being perfect solvers of the parity problem, by not violating the 
constraints derived in the proof, at their various stages of development. 
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